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Abstract — This paper outlines the formulation of a robust fault 
detection and isolation scheme that can precisely detect and 
isolate simultaneous actuator and sensor faults for uncertain 
linear stochastic systems. The given robust fault detection scheme 
based on the discontinuous robust observer approach would be 
able to distinguish between model uncertainties and actuator 
failures and therefore eliminate the problem of false alarms. 
Since the proposed approach involves precise reconstruction of 
sensor faults, it can also be used for sensor fault identification 
and the reconstruction of true outputs from faulty sensor outputs. 
Simulation results presented here validate the effectiveness of the 
robust fault detection and isolation system. 

I. Introduction 

Faults are deviations from the normal behavior of the plant or 
its instrumentation and they can be categorized into: i) additive 
process faults, ii) multiplicative process faults, iii) sensor faults, and 
iv) actuator faults. There exist several fault monitoring procedures 
which can be used to recognize and distinguish different types of 
faults [1], These fault monitoring procedures can be categorized 
into: i) fault detection, ii) fault isolation, and iii) fault identification. 
A survey on design methods for fault detection is given in [2], 
Most of the existing FDI (fault detection and isolation) schemes 
are based on measurement residual generation. Generated residual 
is used to facilitate the decision making procedures involved in FDI. 
The basic difference between most FDI schemes is the underlying 
residual generation methods. Few examples of different FDI schemes 
are the observer based Fault Detection Filters [3]— [5], Kalman filter 
based Proportional Integral Observers [6], Multiple Model Adaptive 
Estimators [7], and system identification methods [8], 

In this manuscript two types of faults are of concern, i.e., actuator 
faults and sensor faults. The FDI scheme considered here is the 
observer based approach. Discontinuous observers such as the sliding 
mode observers have been successfully used in FDI context [9], 
Design of sliding mode observers for detection and reconstruction of 
actuator and sensor faults is presented in [10] and [11], respectively. 
The precise reconstruction of faults proposed in [10] assumes the 
absence of uncertainty. The FDI approaches presented in [11] and 
[12] introduce an approximate fault reconstruction scheme by mini- 
mizing the error between the true fault signal and the reconstructed 
fault signal. In [13], an FDI scheme for a class of nonlinear uncertain 
systems is presented by introducing limitations on the structure of the 
uncertainty. It is important to notice that the precise reconstruction 
of fault is generally not available in the presence of uncertainty. 
While the FDI schemes presented in [10] — [13] involves reconstruction 
of faults, sliding mode observer based FDI approach presented in 
[14] is based on the measurement residual generation. Though the 
FDI scheme presented in [14] assumes precise knowledge of system 
dynamics, a similar FDI scheme which is robust to mismatched 
uncertainties is presented in [15]. The residual generation presented in 
[14] and [15] is based on the sliding mode observer scheme where the 
observer maintain the sliding motion in the presence of mismatched 
uncertainties, but when fault occurs, the sliding motion is broken and 
a residual is generated which contains information regarding the fault. 
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Though there are numerous literatures on fault detection and 
identification, very little work is done on developing robust fault 
detection schemes. Robust actuator/sensor fault detection is a chal- 
lenging problem due to system modeling errors, external disturbances, 
and measurement noise. Most of the current fault detection algorithms 
are model based and they tent to induce false alarms when the plant 
dynamics differ from the assumed model. This manuscript presents 
a robust FDI scheme that can precisely detect and isolate simul- 
taneously occurring actuator faults and sensor faults for uncertain 
linear stochastic systems. The proposed approach is an observer based 
FDI scheme where a discontinuous observer is used for residual 
generation. The main highlights of the proposed FDI scheme are 

• Multiple actuator and sensor faults are considered 

• Proposed technique involves reconstruction of sensor faults and 
therefore this approach can be used for sensor fault identification 
and reconstruction of true outputs 

• There are no constraints on system uncertainties, both matched 
and mismatched uncertainties are considered 

• Present scheme can be easily extended to nonlinear systems by 
considering Lipschitz continuous affine nonlinear terms with 
known Lipschitz constant [13], [15] 

The structure of this paper is as follows. A detailed formulation of 
the observer based FDI scheme is first given. Afterwards, the results 
from numerical simulations and the concluding remarks are given in 
sections III and IV, respectively. 

II. Observer-based Fault Detection Filter 

Let (Q, T, {J~t}t>t 0 , P) denotes a complete filtered probability 
space. Consider an n^-order stochastic system of the following form: 

Xi(f) = AuXi(t) + Ai 2 X 2 (t) + Wi(t) 

X 2 (t) = A 2 iXi(t) + A 22 X 2 (t) + Bu(t) + W 2 (t) 

Yi(t) =CnX 1 (t) + Ci 2 X 2 (t) H-Vr(f) 

Y 2 (t) = C 2 iXr(t) + C 22 X 2 (f) + y e (t) + V 2 (t) 

where Wi(t) and W 2 (f) denote stochastic disturbances and Vi(f) 
and V 2 (t) indicate measurement noises. The state vectors, Xi(t) 
and X 2 (t), are of dimensions, Xi(f) € St"~ r and X 2 (f) € St r , 
respectively. The true state matrices, An E sy^O-OxO-D, a 12 E 

^(n-r)Xr A2i e jjjr: X' (n-r ) ' ^22 £ jjrxr^ and ^ distribu _ 

tion matrix B 6 St rXr , are assumed to be unknown. The desired input 
signal is denoted as u d(t) and u e (t) indicates the error in applied 
control, u(t), due to actuator faults, i.e., 

u(t) = u d (t) + u e (t) (2) 

The stochastic measurement vectors, Yi(t) and Y 2 (i), are of di- 
mensions, Yi(t) € 5R mi and Y 2 (t) € Sf" 2 , respectively. The output 
matrices, Cn, C 12 , C 21 and C 22 , are assumed to be known. The mea- 
surement noise, [V i(t) v^t)] =v(i)er, is assumed to be 
zero-mean Gaussian white noise process, i.e., V(t) ~ A/"(0, RS(t )). 
The vector, y e (f), indicates sensor failures and is modeled as 

y e (f) = f(y e (f), t), Ye (to) = 0 (3) 

where f(-) is an unknown operator. Stochastic external disturbance 
[WT(t) W^(t)] T = W(t) € St" is modeled as a linear system 
driven by a Gaussian white noise process, i.e., 

W(t) =L(W(t),t)+W(t), W(to) = 0 


(4) 



where L(-) is an unknown linear operator and W(t) G is a zero- After appending the sensor error dynamics given in (3), the 
mean Gaussian white noise process, i.e., W(f) ~ A/”(0, Q<5(t)). extended system can be written as 


Assumption 1. Assume the sensor faults are not instantaneous and 
therefore there exist a known conservative upper bound on f(y e (t), t) 
such that 

\f(ye(t),t)\ < a(t), Vf > to 


'Xr(t)‘ 


Amu 

OtAm \ 2 

f3Ami2 

0 ■ 


x 2 „(t) 


1 A 

0,2*77721 

-Am 22 

0 

0 


X 2fl (t) 


0 

0 

Am22 

0 


. y At) . 


0 

0 

0 

^ye. 



where \ ■ denotes the Euclidean norm. 
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The external disturbance, W(t), is mean square bounded [16], 
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where K is a constant matrix whose elements are finite. The assumed 
(known) model of the plant in (1) is given as 


Xmi (() — -(- Aj, 
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Define the model-error vectors T?i(t) G Vt" r and X> 2 (t) G as 

X»i(t) = AAnXi(t) + AAi 2 X 2 (f) + Wr(t) 

X> 2 (t) = AA 21 Xi(t) + A A 22 X 2 (t) + A,Bu d (t) + W 2 (t) 


( 6 ) 


where ATn = An — A mil , AAi 2 = Ai 2 — A mi2 , AA 2 i = 


A 2 i — A m21 , AA 22 = A 22 — A„ 


and A B = B — 


Assumption 2. Given the system parameter uncertainties are 
bounded and the system states are bounded in mean square sense, 
an upper bound on the model error vector T>(t ) can be obtained as 

P (\D(t)\ < fi(t)) = 1, Vt > t 0 

That is, \T>(t)\ is almost surely (a.s.) upper bounded by jX(t) for all 
t > to- 

Now the plant dynamics in (1) can be written in-terms of known 
parameters as 
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where h(-) = f(-) - A ye y e and A ye 
Hurwitz matrix. Let Z(t) = [Xf (t) X 2c< (t) X 2/3 (t) 

now the above extended system can be rewritten as 


™ y!(t)] T , 


Z(t) = FZ(t) + G 3 u d (t) + GiUi(t) + G 2 X» 2 (t) + GsC(t) 

+ G 4 h(-) 


(9) 
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Let H = 


the measurement equations can 
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be rewritten as Y(t) = LtZ(t) + V(t). Now the system in (1) can 
be written as the following dynamically equivalent form 


Z (t) = FZ(t) + G 3 u d (f) + GiUi(f) + G 2 'D 2 {t) 

+ G 3 ^(f) + G 4 h(-) (10) 

Y(t) = HZ(t) + V(t) 


Xi(t) — A mil Xi(f) + A mi2 X 2 (f) -f- T>i(t) 

X 2 (f) — A-m2i Xi(f) + A m22 X 2 (f) + B m Ud(t) + T> 2 (t) + Bu e (t) 

Re-parameterize X 2 (f) as 

X 2 (t)=aX 2o ,(t)+ j 8X 2(S (t) (7) 

where X 2c< (f) G 5R r , X 2/3 (t) G 5R r , a and /3 are user selected scalar 
parameters. Now X 2 (t) can be written as X 2 (t) = q'X 2q (t) + 
/)X 2/3 (f). Select X 2a ,(f) and X 2/3 (f) as 

X 2 „(t) = iA m21 Xi(f) + A m22 X 2a (f) + — X> 2 (t) 

“ 1 ! ° (8) 

X 2j9 (t) = A m22 X 2/3 (f) + -B m u d (t) + —Bu e (t) 


Remark 2. Even though the above representation of the plant is a non- 
minimal realization, the observability of the extended system may be 
obtained by making appropriate changes to the state matrix, F, and 
the corresponding changes to X>i(t), X> 2 (f), and h(-). 

Consider the following partition of Gi as n - r column vectors, 
G 2 as r column vectors, G 3 as r column vectors, and G 4 as m 2 
column vectors as shown below 
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Remark 1. One of the main challenges in the design of observer 
based FDI scheme is the presence of coupled system uncertainties and 
actuator faults [15]. In the presence of coupled system uncertainties 
and actuator faults, it is difficult to design an observer that yields 
measurement residual which is only sensitive to the actuator faults. 
Notice that the re-parametrization of X 2 (t) given in (7) allows 
decoupling of system uncertainties and actuator faults as shown in 
(8). 

Assumption 3. Assume there exists a bounded vector £(t) G 9i r 
such that Bu e (t.) = B m Cit), i.e., 

C(t) = BmBu e (i) and |C(f)| < £(t), Vt > t 0 
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Now the extended system in (9) can be written in summation form 
as 


n — r r 

z(t) =Fz(t) + gi iDi{t ) + g 2j v n - r+j {t) 

7=1 j = l 

r m2 

+ g3k(k(t) + g 4 ihi(-) + G 3 u d (t) 

k=l 1=1 


(11) 



Define G x = [Gi G 2 ] , G 2 = [G 3 G 4 ], and + ft) = 

[£ T (f) h T (-)] T . Now (11) may be rewritten as 

n r+m 2 

Z(t) = FZ(t) + G 3 Ud(t) + y ' GuT>i(t) + Gijtljft) (12) 

i=i j=i 

where Gik and G 2 k are the fc th column vectors of Gi and Gi 
matrices, respectively. Now l = 1, 2, ... ,r + 1 observers of the 
following from are considered 

If t < r 

Z(t ) = FZ\t) + L l [Y(t) - HZ\t)] + G 3 u d (f)+ 

L J (13) 

n Z — 1 r+m 2 

E + E GljVjit) + E GljVj{t) 

i = 1 1=1 l=e+i 

If e = r + 1 

Z*(t) = FZ*(t) + L* [Y(f) - FZ*(t)l + G 3 u d (f)+ 

L J (14) 

71 r + 7712 

E Glil4{t) + E 
i = 1 1=1 

where l/ 6 sjf(n+r+m 2 )xm j s t jj e 0 t, server g a ) n correspond- 
ing to the / !lh observer. The observer inputs are denoted as, 
Hw ••• ^r+m a) (*)] T ^^(*)e» r+m2 . and 

[#*!(*) ••• /4(if = /*'(t)e»”, V€=l, 2, ...,r+l 

Equations (13) and (14) correspond to the typical observer model. 
The observer gain L l and the observer inputs + {f) and u l (t) 
corresponding to the C th observer are selected so that the generated 
residual obtained from observers given in (13) is asymptotically stable 
if there is no fault in the £ th actuator and the residual obtained front 
the observer given in (14) is asymptotically stable despite any actuator 
or sensor fault occurrences. _ 

Define the observer error as Z l (t) = Z(t) — Z e (t). After 
subtracting (14) from (11), the observer error dynamics can be written 
as 

Z{t ) =[F- L l H ] Z l {t ) - L e V(t)+ 

n r+m 2 (15) 

Gu ++) - /4(*)] + [■ ii(t ) - t'j+i 

i = 1 1 = 1 

It is important to note that the solution to the stochastic differential 
equation given in (15) cannot be based on the ordinary mean square 
calculus because the integral involved in the solution depends on 
V(t), which is of unbounded variation, i.e., E [V(f)V T (f + r)] = 
RS(t). For the treatment of this class of problems, the stochastic 
differential equation can be rewritten in Ito form as [18], [19] 

{ 71 
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and therefore the stability of the observer error dynamics given in (16) 
is depicted either as moment stability or stability in probabilistic 
sense. The stability in probabilistic sense is usually known as almost 
sure (a.s.) stability and it is defined as follows [16]: 

Definition 1. The stochastic process Z £ (f) is asymptotically stable 
with probability 1, or almost surely asymptotically stable, if 

P(Z { (t)— >0 as t — > oo) = 1 (17) 

Notice that the almost sure stability of the observer error is 
impossible due to the persistently acting measurement noise S(t). 
Therefore it is desirable for the observer error corresponding to the 
£ <h observer, Z e (t), to have a dynamics that follows 

dZ e m (i) = | [F - L e H ] Z^(t) jdt + L e dB(t) (18) 
Let Z (t) = Z l (t) — Z l m (t), now subtracting (18) from (16) yields 
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Given next is an approach for the selection of the observer gain L l 
and the observer inputs + (t) and u L (t) corresponding to the t h 
observer based on the stochastic Lyapunov approach. 

Since the only information regarding the true observer error is in 

the form of measurement residual, one do not have full access to the 
~ i ~ t ~ i 

signal Z (f), i.e., one only has access to y ( f ) = HZ (f). Based 

~ i 

on (19), dy (f) can be written as 

ay\t) = [f - z\t)+ 

n r+m 2 x 

Y H Gu[Vi(t)-i4(tj\+ E HG 2 j[vAt) - "*+)] \dt 

i= 1 3 = 1 J 

Based on Assumptions 1 and 3, define an upper bound on rj(f) as 
\v(t)\ < v{t), Vf > to 


Theorem 1. Given the Assumptions 1, 2, and 3, the individual fault 
detection filters given in Eq. (13) guarantee almost sure asymptotic 
stability of y (t) if there is no fault occurrence in the £ th actuator 
and the fault detection filter given in Eq. (14) guarantees almost 
sure asymptotic stability of y ( t ) despite any actuator or sensor 
fault occurrences, if the observer gain L e corresponding to the 
observer is selected so that the following matrix Lyapunov inequality 
is satisfied 


[f - Z///j T H T HP e H T H+ 

H T HP e H T H [f - I/#] + Q e < 0 


( 20 ) 


and the observer inputs corresponding to the £ th observer are selected 
as 


where the zero-mean Gaussian white noise V(f) is written as the 
increments of stationary Wiener process with zero-mean and the 
correlation of increments 

F [{S(r) - B(C)} (B(r) - S(C)} T ] = R\t - C| 

Details on stochastic Ito calculus can be_found in [19]. The observer 
error corresponding to the ( th observer, Z t {t), is a stochastic process 


l4(t) = sgn { (y e (t)) T HP e H T HGii } jS(t), 

V* = 1, . . . , n 

Vj(t) = sgn^[y\t)Y HP e H T H G P(t), 

Vj = 1, r + m 2 


( 21 ) 

( 22 ) 



where P l £ K („ +r . +m2 )x(n+r+m 2 ) and Ql € 

^(n+r+m 2 )x(n+r+m 2 ) are p 0S id ve definite symmetric matrices and 
sgn{-} denotes the signum function or the sign function. 

Proof: Construct a Lyapunov function candidate of the form 

V(y\t)) = (y\t)\ HP l H T y\t). Now using the Ito for- 

— £ 

mula [19], dVify (t)) can be calculated as 

dV(y\t)) = l (y\t)Y { [F-L e H] T H T HP e H T H+ 


from the observer given in (14) can be directly used for sensor fault 
detection. That is, if y e (t) = 0, then there is no sensor fault and 
if y ei (t ) 0, then the nonzero y £i ( t ) indicates a fault occurrence 

in the i lh sensor. Moreover, by subtracting y e (t) from the measured 
output yields the true system output. 

III. Numerical Simulations 

Numerical simulation results are presented in this section to 
validate the efficiency of the proposed FDI scheme. Consider a 
stochastic system of the form given in (1) where the true system 
matrices are given as 
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After substituting (20), £V (y e ) can be written as 
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and the assumed system matrices are 
A 

limn — 


Am 2 i — 


0 0 
0 0 


Am 12 — 


1 0 
0 1 


-1.13 0 

0 -1.18 


Am 2 2 — 


B m = 

-0.99 

0 


3.19 

-0.1 

0 

-2.98 


-0.31 

2.01 


The system output 
hx 2 02 x 2 ] , C 12 
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purposes, the external disturbance is mode 
and W 2 (f ) = [W 2 i(t) W 22 (t)] T is given as 
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where the operator £{■} acting on U(x,t) is defined as 

£V(x, t) = bin j^E[dV(X.{t), t)|X(t) = x] (23) 
Substituting (21) and (22) yields 

£V{y e ) < - (y*(t)^ Q e y e (t)+ 

n ( j, 

2 Ej (v*®) HP l H T HGuVi{t)- 


W 2 i{t) = - W 2 i(t) + Wi(f) 

w 22 (t) = - w 22 (t) + m(t) 

where [Wi(f) W2(f)] T = W(f) is zero-mean Gaussian white noise 
process with 

E [w(f)W T (t + t)] = 10" 2 x I 2x2 S(r) 


The measurement noise, V(f) = V(f) £ I 
mean Gaussian white noise process with 


is assumed to be zero- 
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Thus 


Note that T> 1 = 0 and P> 2 {t) is given as 

■D 2 {t) = AA 2 iXi(f) + AA 22 X 2 (f) + ABu d (t) + W 2 (f) 

where the desired control input u d (t) is given in Fig. 1. For the 
re-parametrization of the system states, the constants a and are 
selected as a = /3 = 1. Note that the two possible sensors faults are 


£V(y e (t))<-{y e (t)) T Q e y e (t) (24) 

Therefore the £ = (r + l) th observer in (14) is almost surely asymp- 
totically stable despite the occurrence of any actuator or sensor faults. 
Based on the given proof one could easily make the argument that if 
there is no fault occurrence in the f th actuator, where 1 < £ < r, then 
the £ th observer given in (13) is almost surely asymptotically stable. 
Thus any obser\’ed residual y*(t), will indicate a fault occurrence in 
the £ th actuator. ■ 

Any observed residual in the ( ;|h observer given in (13), where 
1 < £ < r, indicates a fault occurrence in the f lh actuator. Based on 
the observability condition one could easily show that the estimated 
or the observer generated sensor error terms, y e (f), asymptotically 
approaches the true sensor error, y e (f). Therefore y e (t) obtained 



Fig. 1. Desired Control Input 


associated with the fifth and sixth outputs. Two fault scenarios are 
considered here and the details on the fault scenarios are 


1) Fault Scenario I: For the first fault scenario, the faults are 
associated with the second actuator and the fifth output sensor. 
The actuator fault occurs at thirty seconds (sec) and the sensor 
fault occurs at sixty-five sec. Given in Fig. 2 are the u e (t) and 
y e (t) corresponding to the first fault scenario. 

2) Fault Scenario II: For the second fault scenario considered, it 
is assumed that the faults are associated with the first actuator, 
the fifth and the sixth output sensors. The fault associated 
with the first actuator occurs at thirty sec. The fifth and sixth 
output sensor faults occur at sixty-five sec and eighty-five 
sec, respectively. Given in Fig. 3 are the u e (t) and y e (t) 
corresponding to the second fault scenario. 
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Fig. 2. Fault Scenario I 
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Fig. 3. Fault Scenario II 


Note that for the system considered here, there are two actuators 
and therefore three different observers are designed. For both fault 
scenarios considered, A v is selected as A y = 02x2, P e is selected 
as P e = 10 -2 x I, W € {1, 2, 3}. The observer gain is calculated 
as L 1 = L 2 = L 3 = 
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The extended output matrix H and the matrix G can be calculated 
as 
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Details on the results obtained for both fault scenarios are given next. 

1 ) Fault Scenario I: Given in Fig. 4 are the true error vectors, 
2>a(t) = AAaiXr(t) + AA 22 X 2 (t) + ABu d (t) + W 2 (f), C(f) = 
B m 1 Bu e (f), and h(-) = f(-) corresponding to the first fault scenario. 


For simulation purposes the upper bounds on X> 2 (t), £(t), and h(-) 
are selected as 

|D 2l (f)|<10 \T>2 2 (£)| < 20 

ICi(*)l<4 |Ca(*)l < 3 

|/ii(f)|<2 \h 2 (t)\<2 



Fig. 4. Fault Scenario I: Error Vectors 




(b) Estimated Sensor Faults 


Fig. 5. Fault Scenario I: Observer Residual and Estimated Sensor Faults 


Given in Fig. 5 are the generated residual and the estimated sensor 
errors corresponding to the first fault scenario. Figure 5(a) contains 
the measurement residual generated for observer one and two. The 
first two kinks in the residual are due to the start of input application 
that occurs around five seconds and the leveling-off the input to its 
steady state value around twenty seconds. Notice the jump in observer 
two residual around thirty seconds due to the fault occurrence in 
the second actuator. Figure 5(b) contains the estimated sensor errors 
obtained from the third observer. Note that the estimated sensor error 
is similar to the true sensor error given in Fig. 2(b). 

2) Fault Scenario II: Given in Fig. 6 are the error vectors, 
X> 2 (t), £(f), and h(') corresponding to the second fault scenario. 
The upper bounds on error vectors used here are the same upper 
bounds used for the first fault scenario. 

Given in Fig. 7 are the generated residual and the estimated sensor 
errors corresponding to the second fault scenario. Figure 7(a) contains 
the measurement residual generated for observer one and two. Notice 
the sudden increase in observer one residual around thirty seconds 
due to the fault occurrence in the first actuator. Figure 7(b) contains 
the estimated sensor errors obtained from the third observer. Note 
that the estimated sensor errors are similar to the true sensor errors 
given in Fig. 3(b). 



Fig. 6. Fault Scenario II: Error Vectors 




(b) Estimated Sensor Faults 


Fig. 7. Fault Scenario II: Observer Residual and Estimated Sensor Faults 


IV. Conclusion 

Robust actuator/sensor fault detection is a challenging problem due 
to the effects of modeling errors, system process noise, and mea- 
surement noise. This manuscript outlines the formulation of a robust 
fault detection and isolation scheme that can precisely detect and 
isolate simultaneously occurring actuator faults and sensor faults for 
uncertain linear stochastic systems. The given robust fault detection 
scheme would be able to distinguish between model uncertainties 
and actuator failure and therefore eliminate the problem of false 
alarms. The presented approach involves precise reconstruction of 
sensor faults and therefore this approach can be used for sensor fault 
identification and the reconstruction of true outputs from faulty sensor 
outputs. The proposed approach is an observer based fault detection 
and isolation scheme where a discontinuous observer is used for 
residual generation. 

The proposed approach assume conservative upper bounds on the 
system uncertainties, the actuator faults, and the sensor fault rates. 
A bank of discontinuous observers is designed for fault detection 
and isolation scheme where the number of observers is based on 
the number of actuators. The observer gain and the discontinuous 
observer inputs are selected so that the observed residual is almost 
surely asymptotically stable if there is no actuator fault occurrence. 
As a result, any observed residual would indicate a fault occurrence 
in the corresponding actuator. In addition to the bank of observers 
designed for actuator fault detection, a robust discontinuous observer 
is designed so that the estimated or the observer generated sensor 
error terms asymptotically approaches the true sensor error. There- 
fore, the sensor error estimates obtained from the robust observer 


can be directly used for sensor fault detection, isolation and identi- 
fication. Moreover, by subtracting the estimated sensor errors from 
the measured outputs, true system outputs can be generated. The 
simulation results reveal clear indication of actuator faults despite the 
presence of matched system uncertainties and external disturbances. 
Moreover, the estimated sensor errors are identical to the true sensor 
error regardless of the measurement noise present. 
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